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We review here the development of the general formalism for the study of fermion propagation 
in the presence of stochastic media. This formalism allows the systematic derivation of evolution 
equations for averaged quantities as survival probabilities and higher order distribution moments. 
The formalism equally applies to any finite dimensional Schrodinger equation in the presence of 



a stochastic external field. New integrodifferential equations valid for finite correlated processes 
are obtained for the first time. For the particular case of exponentially correlated processes 
a second order ordinary equation is obtained. As a consequence, the Redfield equation valid 
for Gaussian delta-correlated noise is rederived in a simple way: it has been obtained directly 
and as the zero-order term of an asymptotic expansion in the inverse of the correlation length. 
The formalism, together with the quantum correlation theorem is applied to the computation of 
higher moments and correlation functions. It is shown that equal and unequal time correlators 
follow similar differential equations. 

X 

1. INTRODUCTION 

We review here the general formalism for the study of neutrino propagation in the presence 
of stochastic media [ |l| . This formalism will allow the systematic derivation of evolution equa- 
tions for averaged quantities as survival probabilities and higher order moments of them. New 
integrodifferential equations valid for finite correlation processes are obtained for the first time. 
For exponentially correlated processes a second order ordinary equation is obtained as a conse- 
quence. The Redfield equation valid for Gaussian delta-correlated noise is rederived in a simple 
way as a particular case: it can be obtained directly or as the zero-order term of an asymptotic 
expansion in the inverse of the correlation length. In the context of neutrino oscillations, it is 
shown that the presence of matter noise induces the appearance of a effective complex part in the 
matter density. Finally the formalism, together with the quantum correlation theorem, will be 
applied to the problem of the computation of distribution higher moments and non-equal-time 
correlation functions. 

The formalism may be applied to any quantum system governed by similar Schrodinger equa- 
tions: equations where a stochastic function appears multiplicatively in some part of the Hamil- 
tonian. For simplicity of notation, in this work we will always deal with Hamiltonians which are 
finite dimensional operators in some Hilbert space, the equations will be obviously of the same 
type for general systems described by infinite-dimensional Hamiltonians. 
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2. THE GENERAL CASE 

Let us consider a general system whose evolution is described by the following linear random 
Schrodinger equation: 

id t X = p(t)L(t)X, X(0) = X o , (1) 

where X is a vector (i.e., a multiflavor wave function) or a matrix (i.e., a density matrix) of 
arbitrary dimension. Equation ([!]) is defined in an interaction representation, any additive non- 
stochastic term has been solved for and absorbed in the definitions of X(t) and L(t). L(t) is a 
general linear operator. It can be an ordinary matrix if X represents a wave function. If X is a 
density matrix, L(t) is a commutator: 

L A (t)X = [A(t),X]. (2) 

Any other linear operator is admissible: obviously it is always possible to define an enlarged 
vector space where the action of L is represented by a matrix. 

We will assume that p(t) in Eq.© is a scalar Gaussian process completely determined by its 
first two moments, which, without loss of generality can be taken as 

(p(t))=0, (p(tW)) = f(t,t'). (3) 

In other terms, p(t) is characterized completely by the measure 
1 f°° 

[dp]exp--/ P (t)r\t,t')p{t')dtdt'. (4) 

An important particular case is when the process is 5-correlated, the correlation function is 
then of the form 

f(t,t') = n 2 5(\t-t'\). (5) 

A convenient way of parametrizing a correlation function with a finite correlation length is to 
use an exponential function: 

f{t,t') =0 2 eexp(-e | t - t' |) . (6) 

In this case the correlation length is defined as r = 1/e. The expression (ph is reobtained 
letting t — > or c — ► oo in Eq.(|6|). 

X(t), the solution to Eq.(g), is a stochastic function. The objective of this work is to obtain 
equations for its ensemble average and higher moments in a systematic way . Let us remark 
that for us Eq.(|l|) is purely phenomenological, we assume that it is the result of a more complete 
microscopic analysis which can account for the randomness of p(t) (see for example, Refs.[ 

i 10,11 13)- 



To obtain a differential equation for the average of X we make use of the following well known 
property: For any Gaussian process p(t) characterized by a <5-correlation function (||) and any 



functional F[p] we have the following relation [ 12 



For a Gaussian process with an arbitrary correlation function [Eq. (Q)] we have instead the 
general relation: 

(F[p]p(t)) =Jdr (p(t)p(r)) . (8) 
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Let us consider now the evolution operator U for a particular realization of Eq.(Q). By 
definition, 

X(t) = U(t,t )X . (9) 

The operator U(t, to) is a functional of p, it has the following formal expression in terms of a 
time ordered exponential: 



U(t, to) = Texp 



to 



p(t)L(t)cIt 



(10) 



The functional derivative of U with respect to p{t) can be computed by direct methods. By 
differentiating term by term the series expansion for Eq.(|i~0|), the result is 

6 JL^1 = -iL(r)U(r,t Q ), t<r<t . (11) 

In order to obtain a differential equation for {X) we observe that 

i(8 t X) =id t (X) =L(t)(p(t)X). (12) 

Combining together Eqs. (P1)>(PD>(PD we easily obtain the result that the evolution equation for 
the ensemble average is in the general case an integrodifferential equation given by 

id t {X(t)) = -i tdt 1 f(t,t')L(t)L(t'){X(t')), (X(0)) =X Q . (13) 



o 

This equation, which is exact and of very general validity, is the equation which we were looking 
for and one of the main results of the present work. Note that previously integrodifferential 
equations have been obtained, valid for particular cases or in particular limits, using heuristic 
ad hoc arguments (for example, in Ref. [ |l(J). The derivation of Eq.([l3|) which has been done 
here is the rigorous justification for such approaches. 

In some notable cases Eq.(13) can be reduced to an ordinary differential equation. In the next 
sections we will see how it can be reduced to, respectively, first and second order equations for 
the correlation functions given by Eqs. @ and (|). 

3. THE PARTICULAR <5-CORRELATED CASE 

For the particular case where the correlation function is of exponential type, a second order 
ordinary differential equation can be obtained as we will see below. On the other hand, for the 
simpler <5-correlated case the evolution equation (13) becomes the ordinary differential equation 

id t (X(t)) = -iQ 2 L 2 (t)(X(t)). (14) 

Taking L as a commutator, this last equation coincides with the Redfield equation derived by 
Ref. [ 9. Note that the effective " Hamiltonian" appearing in the second part of Eq. (|i~4|) is not 
Hermitic anymore (an example of a fluctuation-dissipation effect). 

In cases of interest for the neutrino oscillation problem, the original equation for the density 
matrix is of the slightly simpler form 

id t X = [Ho(t)+p{t)g(t)H u X], X(0)=X , (15) 



where we have written the commutator explicitly. p(t) is a stochastic function as before and 
g(t) an arbitrary scalar function. Hq,H\ are Hamiltonian matrices, the former contains the 
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average part of p(t): Ho = H Q (t) + po(t)Hi, with po = (p). The latter is assumed to be time 
independent. 

For the problem described by Eq. (|l5|) the corresponding Redfield equation for the averaged 
density matrix is of the form 

id t (x) = [H (t),{x)]-in 2 g 2 (t)[H 1 ,[H 1 ,(x)]], (is) 

= H (X)-(X)H+ + 2in 2 g 2 (t)H 1 (X)H 1 , (17) 
where in the last line the following effective Hamiltonians were defined: 

H± = H ±ig 2 (t)H 2 . 

The solution of what is called the "coherent" part of Eq. (|T^ ) (two first terms of the Hamil- 
tonian) is accomplished by defining the average evolution operator: 

(C7±) = TexpjdrH^r), (U~) = . (18) 

The coherent part of the density matrix is then 

Wcoh = {U-)X {U+)K (19) 
Defining a new "coherent" interaction representation by the relations 

H L = (U-y'^iU-), H R = (U-)H 1 (U-)-\ (X) I = (U-)(X) coh {U-)-\ (20) 

we arrive at the equation we were looking for, the resolution of the original equation is equivalent 
to the resolution of the following one: 

id t (X)j = 2m 2 g 2 (t)H L (X)jH R . (21) 

There are some important particular cases where Eq. ([i~4|) can be solved or considerably 
simplified by taking into account the algebraic properties of a specific L (in what follows k(t) is 
always a scalar function). 

(a) Let us assume that L is such that L 2 (t) = k(t)L(t). This case appears in the computation 
of the average wave function with matter density noise. Equation ([14] ) reduces to 

id t (X) = -in 2 k(t)L(t)(X). 

The averaged equation is similar to the original one, the nonrandom part of the density is 
"renormalized" acquiring an imaginary term 

P — > Po — iQ 2 k. 

This is the density which will appear in the coherent effective Hamiltonians Hq . 

(b) The case where L 2 (t) = k{t)I with I the identity matrix appears in the computation of 
the averaged neutrino wave function under noisy magnetic spin-flavor precession. The resulting 
equation can trivially be integrated (to be compared with the previous case): 



(X(t)) = exp 



-n 2 f dTk(r) 
Jo 



The average wave oscillation is damped by a factor equivalent to the one first calculated by 
Nicolaidis [ A similar damping also appears when computing the average density matrix 
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from Eq. (|2l|). From these differences of behavior with respect to case (a) it is expected that the 
presence of magnetic field noise can have some influence even if applied far from any resonance 
region. 

(c) The case where L (t) = —k(t)L 2 (t) appears in the computation of the average density 
matrix with matter or magnetic noise. We can obtain in this case the "conservation law" 

l-fl 2 k{t)] L 2 (t)d t (X{t)) =0. 

L 2 (t) is not invertible because the operator L(t) has a zero eigenvalue. The previous expression 
has proved to be of practical importance in some concrete numerical applications [ 13, 14, |l5| . 

The presence of zero eigenvalues distinguishes cases (a) and (b) from (c); it can have conse- 
quences in the long term behavior of the respective ensemble averages as is shown elsewhere [ 



4. AN ASYMPTOTIC EXPANSION FOR EXPONENTIALLY CORRELATED 
SYSTEMS 

We will see now how Eq. (|D can be obtained as a limiting case when the correlation length 
tends to zero. For this purpose we use an exponential correlation function as Eq. (^), the 
integrodifferential evolution equation becomes in this case 



idt{X(t)) 



-m 2 eexp(-et) /* dt' exp(eT)L(t)L(t')(X(t')). 
Jo 



(22) 



Let us compute the asymptotic expansion of the second term of Eq. valid for large e; the 
following expansion is valid for any function g(t): 



He) 



eexp(— et) / dr exp(er)g(r) ~ g(t) 
Jo 



g\t) , g"(t) 



+ 



+ 



Inserting this expression into Eq. 



idt(X) 



SI 



in 2 L 2 (t)(X) + i—L(t)d t (L(t)(X)) + o[- 



we obtain the following expansion in powers of e: 



1 



(23) 



(24) 



To leading order in 1/e, we recover the expression corresponding to the <5-correlated case. At 
next-to-leading order we get finite-correlation correction terms 

n 2 „. .<>- 



idt(X) = 
or, equivalently, 



itt 2 L 2 (t)(X) + i—L(t)L'(t){X) + i—L 2 (t)d t {X) 



1 - — L 2 (t) 



dt(X) 



n 2 

-n 2 L 2 (t) + —L{t)L'(t) 



(X). 



(25) 



(26) 



Finally, We get the following differential equation valid to order 1/e, making the assumption 
that the operator which multiplies the left term is invertible: 



dt(X) 



-n 2 L 2 (t) + —L(t)L'(t) + — L\t)) (X). 
e e 



(27) 



This equation can be used for finite, but relatively large, correlation lengths. We see that, up 
to this degree of approximation, not only the ratio level of noise to correlation length (J7 2 /e) is 
important. We have different regimes according to the value of Q 2 . The first term will be more 
important for a low noise amplitude (ft 2 « 1). For strong noise (Q 2 » 1) the second term, 
proportional in this case to L 4 , will dominate. 
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5. AN EXACT DIFFERENTIAL EQUATION FOR EXPONENTIALLY CORRE- 
LATED SYSTEMS 

In contrast with the approximate approach used in the previous section, we can actually derive 
a simple, ordinary second order differential equation for the case of exponential correlation. Let 
us assume that the original equation is of the same decomposable type as the one appearing in 
Eq. ( |l5| ) but let us include other cases using the general notation 

id t X = (Lo(t) + p(%(t)£i)*, X(0)=X o , (28) 

with Lo,Li general linear operators as before, the latter time independent. In this case Eq. (|l3|) 
is of the form 

id t (X) = L (t)(X) - iQ 2 eLle- et g(t) f dt'g(t')e et ' (X(t')). (29) 

Jo 

Differentiating the equation once and performing some simple algebra we obtain the following 
ordinary second order differential equation: 

[d t -X(t)][id t -L (t)}(X) = -in 2 eg 2 (t)Lj(X), 

i(d t (X)) = L(0){X) o , (X) Q = X , (30) 

with 

\(t) = -e + g'(t)/g(t). 

Let us remark that this equation is exact and probably the most important result of this work 
from a practical point of view. 

6. HIGHER ORDER MOMENTS AND THE QUANTUM REGRESSION THE- 
OREM 

Second order distribution moments, expressions of the type (X^Xj), or in general, moments 
of any order, can also be computed using equations similar to Eq. (|13|). The straightforward 
procedure is to define products X^...^ = XiXj ■ ■ ■ Xk and write differential equations for them 
using the constitutive equations for each of the Xi . The resulting equations are of the same type 
as Eq. (jl|). The similarity is obvious when one adopts a tensorial notation and defines products 
of the form X^ = X ® • • • (g) X. The task of obtaining evolution equations is especially simple 
within this notation. 

Correlators of quantities at different times also appear in the computation of averages of 
quantities of physical interest, i.e., expected signal rates. In the most simple case expressions 
of the type (Xi(t + T)Xj(t)). We will shortly show that the knowledge of equal-time moments 
and the application of the quantum correlation theorem is sufficient for the computation of this 
kind of correlators. 

Our objective now is the computation of the correlation function appearing in two steps. First 
we define the generalized density matrix X^ as the tensorial product of usual density matrices 
at two different times and energies: 

xW(E 1 ,E 2] t 1 ,t 2 ) = X(E 1 ,t 1 )®X(E 2 ,t 2 ). (31) 

(2) 

The average of the element {X^'^) = {Xi±X±i} is evidently the probability correlation function 
we are looking for. 
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The differential equation for the equal time function X^ {E\, E% t, t) is obtained from the 
individual evolution equations for the matrices X\p = X(E\p) [indices (1,2) label expressions 
where E\,E<i appear respectively] 

d t X® = H X X® +xWH 2 = pLX®. (32) 



Equation (|32| ) is a random linear differential equation, linear in the stochastic variable p(t). 
Applying the formalism developed in the previous section, we can immediately write the equation 
for the ensemble average (X^) [Eqs. ([l3]-[l4])]. Once we know the equal time correlator, we 
obtain the expression for any other pair (t, t') using the quantum regression theorem [ 15] which 
reads as follows. For any vector Markov process Y, if the ensemble average of Y fulfills a 
Schrodinger-like equation of the type 

d t (Y(t)) = G(t)(Y(t)), (33) 

with G(t) an arbitrary matrix, then the second order correlations will obey the following equa- 
tion: 

d T {Y i {t + T)Y l {t)) = Y J G ij {r){Y j {t + T)Y l {t)). (34) 

3 

Note that the quantum regression theorem is in principle not applicable to systems described 
by the integral equation Eq. (^). Nevertheless it is applicable to problems where the integral 
equation can be reduced to an ordinary differential equation such as Eq. (|30|) . Such a second 
order equation can easily be expressed as a first order one by defining the auxiliary pair process 
Y = (X, d t X). 

For the simplest case where p is a 5 correlated process the equation for non-equal time corre- 
lators is explicitly given by 

id T (X^{t + T,t)) = -ifl 2 L 2 (T)(xW(t + T,t)). (35) 

The initial condition, X^(t,t), is the equal time correlator previously obtained. For expo- 
nentially correlated problems second order equations similar to Eq. ( |30| ) can immediately be 
obtained. 



7. CONCLUSIONS AND FINAL REMARKS 

In conclusion, in this work we have developed the general formalism for the study of neutrino 
propagation in stochastic media. This formalism has allowed the systematic derivation of evo- 
lution equations for averaged quantities as survival probabilities and higher order moments of 
them. New integrodifferential equations valid for finite correlation processes have been obtained 
for the first time. For exponentially correlated processes a second order ordinary equation is 
obtained as a consequence. The Redfield equation valid for Gaussian <5-correlated noise is red- 
erived in a simple way: it has been obtained directly and as the zero-order term of an asymptotic 
expansion in the inverse of the correlation length. 

The formalism can be generalized in an obvious way to obtain the ensemble average of equa- 
tions of a slightly more general type than Eq. ([!]), equations of the form 

id t X = ^pi^Li^X, X(0) = X . (36) 

i 

These equations appear, for example, in the chaotic neutrino magnetic precession. The result 
for the case where the pi are <5-correlated in time but mutually uncorrelated is simply 

id t (X(t)) = -iJ2^L 2 (t)(X(t)). (37) 
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The generalization to continuous Schrodinger equations in the presence of random potentials or 
random external forces is also obvious if we leave aside the mathematical differences coming from 
the appearance of infinite-dimensional Hilbert spaces. The formalism is of general application to 
any quantum system governed by similar Schrodinger equations: equations where a stochastic 
function appears multiplicatively in some term of the Hamiltonian. 
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